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SOME GENERALIZED POWER SERIES INVERSIONS 

R. H. Estes 
E. R. Lancaster 

ABSTRACT 

A function G(£, v, t) = 0 is expanded as a tri-variate power 
series and a general inversion algorithm derived. Various im- 
portant special cases are examined, including a two dimensional 
generalization of the Schlbmilch- Cesar o recursion formula for the 
high order derivatives of a function of a function. 
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SOME GENERALIZED POWER SERIES INVERSION 


INTRODUCTION 

The use of recursive power series on digital computers as a numerical tech- 
nique is well known in the literature. Moreover, with the development of com- 
pilers for the formal manipulation of symbolic expressions the method of re- 
cursive power series has also become an important tool in analytical studies. 
However, most efforts in these fields have been directed toward efficient algo- 
rithms for various highly specialized problems. It is the purpose of this paper 
to present a very powerful, general algorithm which itself may be specialized 
to solve a variety of specific problems. 

A GENERAL TWO DIMENSIONAL INVERSION 
Consider a function of three variables 

G(<f, 7J, T) = 0 

g (£o-V t o) = 0 

which may be represented by a power series in a region about f 0 , t) 0 , t 0 . 
let 


( 1 ) 

Then 


y v - V 0 t=r-r 0 

so that 

x o = y„ = l o = 0 


i 


and 


G( x +£ 0 , y + V t +r 0 ) = F(x, y, t) = 0. 


( 2 ) 


Then 


0 = FXx, y, t) 


00 00 00 

a. . x 1 y* t k 

i , j , k J 

i = 0 j = 0 k=0 




i=l j = 0 k=j 


( 3 ) 


where 


i ! j ! k! a 


<j i+j + k 


ijk 


<?x‘ dy 1 dt 


k F(x, y, t) 


a ooo 


0 . 


Put 


W. 

I 


S x 1 



W. 


ik 


t k 


v , = y* = £ v lk t» . 

k= i 


( 4 ) 


Writing 


W : 


"i-i 


V; = W, 
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and forming the Cauchy product yields the recurrence formulas 


i-i+l 


W. . 

i j 


21 

k = 1 


V.. 

1 J 


T 

k=l 


V V 

V lk v i-l,j-k 


( 5 ) 


Note that 


W 0j V Qj “ 8 0j ; W jk 0 , j > ^ 

where S.. is the well-known Kronecker delta. Substituting Equations (4) into (3) 
we obtain 


0 s F(x, y, t) 


00 00 00 



i = 0 j=0 k-i + j 


( 6 ) 


where 


'ijk 


k-i-j 

L 

m- 0 


a. . G - , * 


G 


ijk 


£ 


w. . V. 

i , k~ n ) n 


( 7 ) 


Observe that 


^Ojk ^jk ®iOk ®ik ®00k ^Ok ' 
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Equating to zero the coefficients of powers of t in Equation (6) yields 


i i-k 


^ ^ j » ’ i - 0, 1, 2, 


k*0 j=0 


Using Equations (7) we have for i = 0, 1, 2, 


i i~k i-k- j 


**kj m ^k , j , i~n 


k=0 j = 0 m s 0 


i i~ ^ i~ k~ j i“m~k 


jm^k , i-m-n 


k=0 J = 0 tn=G n=j 


Rearranging (9) gives, for i = 0, 1, 2, 


i i 

( a 010^1i + a i00*li) ~ a 00i + y ’ [ a 01m^l,i-ni + a iQm^l,i-iii] 


i i~ J 

7 "! H [ a oi« v i.i-« +a io» w i.i^] 


j = 2 m=0 


i“ 2 i-m-1 


rnrn^l.n 


m=0 n=l 


i“l i-j-l i~m~ 1 


, i“nrn ^ j, n 


j = 2 m*0 n* j 
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+ 


i-l i-j-1 i-m-j 


LLL 

i = 2 m- 0 n- 1 


j lm 


W 

j f i“nr*n 


v, 


♦ n 


i~2 i-k i ~k“j 

•LLL 

k“2 j * 2 m-0 



n=j 


*ki 


W V 

i k, i-nrn j,n 


( 10 ) 


where it is understood that a sum is identically zero if the upper summation 
index is less than the lower index. Hence if x or y is known as a power series 
in t, Equation (10) yields y or x respectively as a power series in t. 

If we are given two functions in the form 


0 s F(x. y- 1 ) 


00 00 00 



i=0 j=0 k=0 


00 00 00 

v — » 


o ■ HCx,y,t) = 2^Z^Z. bijkXiyJtk 

i = 0 j = 0 k*0 


( 11 ) 


where 


a oio^ioo " a ioo^oio ^ ® 


a 000 ^000 ® 


there will be two equations of the form (10) which may be solved simultaneously 
for Vjj and W u , thereby giving 

x = 

i = ! 
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y 


(12) 


oo 

j=» 

The algorithm for computing the W, . , V, . when given a. jk , b. jk is as follows: 
For i = 1, 2, •••, j - 2 compute tj and V J _ i { from Equations (5). Then 
solve the two equations of the form (10) for W t . and . 

To indicate the flexible nature of Equation (2), let 

F(x, y, t) = g(x, y, t) - f(t) 3 0 


where 


oo oo oo 




i=0 j-0 k*0 


f(t> = 2Z fit ‘ 

i s 0 


and where F(x, y , t ) is given as an expansion by Equation (3) and x and y are 
given by Equations (4). Then suppose we are given g. jk , W, ^ and for i, j, k 
= 0, 1, 2, ••• and we desire f.. We have 

8 i;0 ®ij0 

®00k ' ®00k ~ ^k 

a ijk = «ijk • ? 0 
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and Equation (10) becomes, for i = 0, 1, 2, 


* « t 


i-1 

{ i ~ gpoi * 8oiO V ii * glOQ W li + } 

m= 1 


i i“ j 


i“ 2 i-nrl 




j = 2 m=0 


m=0 n= 1 


‘ I i“i-l i-n 


i~ 1 i ~ j ~ 1 i “mr j 


21 21 22 *».•!. i-rs*. * 22 22 22 


j - 2 m= 0 


n= j 


p2 m=0 n= 1 


+ 



k=2 


i”k i“k-j 



j x 2 m^O 


i" m~k 



s k j m 


w. . v. . 

K i i iit^ n ) n 


TWO-DIMENSIONAL SCHLOMILCH-CESARO FORMULA 

An important special case of our general result is the two dimensional 
analogue of the Schlbmilch- Cesaro recursion [1] which gives the total derivative 
of a function of two variables which in turn are functions of a single variable. 
Szebehely [4] has generalized the Faa di Bruno formula [2] to this problem but 
his result is somewhat cumbersome for numerical application. To obtain our 
special case let Equation (2) take the separable form 


F(x, y- t) - h(x, y) - f(t) e 0 

where x and y are functions of t and where 

x o = y 0 = t 0 = 0 • 


( 13 ) 
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Put 


h(x. y) 



Then by Equation (3) 


a i jk 


f(t) 

II 

f. t 1 . 

i 

r o 

for 

i, j, or k / 0 

A ii 

for 

o 

II 

- K 

for 

O 

II 

II 

A - f 

^ ft 00 0 

for 

i = j = k = 0 


(14) 


(15) 


Equation (9) then reduces to, for i = 1, 2, • • • 


1 

k=l 


i-1 i-k i 
k- 1 j=l m-k+j 


(16) 


where our previous summation convention has been used and 


“000 


= \o - f o = 0 
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Then given A^, W tj and V ljt we can calculate f. from Equations (16) and 
(5). 

LOWER DIMENSIONAL SPECIALIZATIONS 

If we let y be identically zero in Equation (2) and sin plify Equation (10) we ob- 
tain a special case of our more general result which itself has many important 
and interesting special cases. Equation (2) becomes 

g(x, t) = 0 

*o = % = 0 (17) 


and y = 0 implies 

v 00 = 1 , V ij = 0 , i, j t o 

Equation (10) reduces to, for i = 1, 2, • • • 




a 0i " 


1 

L 

j = 2 


a ioWji 


££ 

j=l 


V w i,i- k 


(18) 


k= 1 


where a. . replaces a.„. and a rn =0. 

1J r iOj 00 

Equations (17) and (18) give a general result which contains the following special 
cases of interest; the expression of the solution x of the equation f(x) = h(t)as 
a power series in t as derived by Thacher [5] , the recursive inversion of power 
series, and a recursive Schlomilch-Cesaro algorithm to compute the derivatives 
of a single-variate function of a function. 
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To obtain Thacher’s result we assume Equation (17) is separable in the form 


g(x, t) = f(x) - h(t) = 0 


(19) 


where x is a function of t . We then have 


l io Ijli dx j 


a oj 


/ 1 \ d> h 

Vj'.i dt i 


= i 'll g 

'i j Vi! j ! / At j 


<?x‘«?t J 


= 0 , i, j / 0 


Equation (18) becomes 


no 


1 

j = 2 


( 20 ) 


Denoting 


CD 

f(x) = £\ xi 
i=0 


h(t) = 

i = 0 


( 21 ) 


we have 


f j = a io - h i = " a oj 
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i - 2, 3. • • • 


( 22 ) 


f l*U 


- h. 


f .*n 


- h. 


1 

L 

j = 2 




Thus, when f and h are given as power series and the root x = 0 is known at 
t = 0, the solution of Equation (19) is given as 


x - 



For the further specialization 


h(t) - t 


or 


h. = S. 


i 1 


we have the well-known power series inversion problem whereby given the series 


00 



i=l 


we desire 


x - 



*n 


t ' 1 . 
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The need to compute high order derivatives of a function of a function 


f(x) - f[x(t)] - h(t) 


(23) 


occurs frequently. Let 


f(x) = £ f i xi 
i=l 


i=l 


h(t) = 

: = i 


(24) 


Then given f . and W u , we desire h . . An algorithm presented by both Schlomilch 
and Cesaro [3] gives 


1 

r — "l 

** = 

j=i 


where P. . is a polynomial in W lk , k = 1, • • • , i , established by the formula 


P,.U : TtPi.k • k= 1.2. -.i*l 


where 


P 


i + 1,0 


0 , 


p 


i+l ,k 


0 . 


k > i+l . 
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However, Equation (22) may be written 


i 

h, = £ f ' W i‘ <25) 

where the W j . are determined recursively in terms of the W n from Equation (5) 
and no further differentiation is necessary. 
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